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Abstract 

The ground state and magnetization process of the mixed spin-( 1,1/2) Ising diamond chain is exactly solved by employing the 
generalized decoration-iteration mapping transformation and the transfer-matrix method. The decoration-iteration transformation 
is first used in order to establish a rigorous mapping equivalence with the corresponding spin-1 Blume-Emery-Griffiths chain in a 
non-zero magnetic field, which is subsequently exactly treated within the framework of the transfer-matrix technique. It is shown 
that the ground-state phase diagram includes just four different ground states and the low-temperature magnetization curve may 
exhibit an intermediate plateau precisely at one half of the saturation magnetization. Our rigorous results disprove recent Monte 
Carlo simulations of Zihua Xin et al. [Z. Xin, S. Chen, C. Zhang, J. Magn. Magn. Mater. 324 (2012) 3704], which imply an 
existence of the other magnetization plateaus at 0.283 and 0.426 of the saturation magnetization. 
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magnetization plateau in the low-temperature magnetization 
process To provide a comprehensive descrip- 



1. Introduction 

During the last few years, a considerable research interest 
has been devoted to the frustrated magnetism of diamond spin 
chains that was originally initiated by the effort to clarify sev- 
eral unusual magnetic features of the natural mineral azurite 
Qi3(C03)2(OH)2 like for instance a presence of the one-third 
zation plateau i 
flilii. 

tion of the overall magnetic behaviour of the azurite, which is 
basically affected by a mutual interplay between the geometric 
spin frustration and quantum fluctuations, one has to employ a 
rather sophisticated combination of the first-principle density- 
functional calculations supplemented with the extensive numer- 
ical calculations 

aa. 

However, it is worthy of notice that 
some important vestiges of this intriguing magnetic behaviour 
can be traced back also from the relevant behaviour of the ex- 
actly tractable spin- 1/2 Ising diamond chain 0]. 

On the other hand, the exactly solvable mixed-spin Ising di- 
amond chains have received much less attention so far. To the 
best of our knowledge, the mixed spin-(l/2,l) Ising chains with 



the spin- 1/2 nodal atoms and the spin-1 interstitial (decorating) 
atoms have been just marginally investigated as special liming 
cases of the exactly solved mixed-spin Ising-Heisenberg dia- 
mond chains J3,SSl- Quite recently, Zihua Xin et al. [ 10] have 
studied another version of the mixed spin-( 1,1/2) Ising diamond 
chain with the spin-1 nodal atoms and the spin- 1/2 interstitial 
atoms within the framework of numerical Monte Carlo simula- 
tions. The main purpose of this work is to provide the exact an- 
alytical solution for this mixed-spin Ising diamond chain, which 
will convincingly contradict a presence of the striking magneti- 
zation plateaus at 0.283 and 0.426 of the saturation magnetiza- 
tion theoretically predicted in Ref. ifioll . 

This paper is organized as follows. The model and basic steps 
of the exact method will be clarified in Sec. [2] The most inter- 
esting results for the ground-state phase diagrams and magneti- 
zation process of the symmetric as well as asymmetric mixed- 
spin Ising diamond chain will be widely discussed in Sec. [3] 
Finally, our paper will end up with several concluding remarks 
and future outlooks mentioned in Sec. |4] 



2. Model and its exact solution 
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Let us begin by considering the mixed spin-( 1,1/2) Ising di- 
amond chain in a presence of the external magnetic field. The 
magnetic structure of the investigated model system is schemat- 
ically illustrated in Fig. Q] together with its primitive unit cell. 
As one can see, the primitive unit cell in a shape of diamond 
spin cluster involves two nodal Ising spins Sk and Sk+i along 
with two interstitial Ising spins yiztj and 11^,2- The total Hamil- 
tonian for the mixed spin-( 1,1/2) Ising diamond chain in a pres- 
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Figure 1: A fragment from the mixed-spin Ising diamond chain. The nodal (Sic, 
Sk+i) and interstitial (/i^i, fik.l) Ising spins belonging to the kth primitive unit 
cell are marked. The Ising interactions J[ and 73 along sides of the diamond 
unit cell are equal (different) in the symmetric (asymmetric) diamond chain. 



ence of the external magnetic field reads 
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where the nodal Ising spins Sk - ±1,0, the interstitial Ising 
spins /j.k. a - ±1/2 (a = 1,2) and the periodic boundary condi- 
tions yUo,i = Hn,\,Ho,2 — f*N,2 are implied for convenience. The 
interaction constants Ji and ^3 label the nearest-neighbour in- 
teractions between the nodal and interstitial Ising spins along 
sides of the primitive diamond unit cell, while the interac- 
tion term Jo accounts for the diagonal interaction between 
the nearest-neighbour interstitial spins from the same primitive 
cell. Finally, the Zeeman's term h determines the magnetostatic 
energy of the nodal and interstitial Ising spins in the external 
magnetic field. 

For further manipulations, it is quite advisable to rewrite the 
total Hamiltonian (Q3 as a sum over cell Hamiltonians 



(2) 



k=\ 



whereas the cell Hamiltonian < H k involves all the interaction 
terms of the kth diamond unit cell 

<H k = Ji(jJk,iSk + Hk,iSk+\) + JiUKiHk,! + Ji(Hk,\Sk+\ +Hk,iSk) 



h 

-Kix k ,\ +lXk,i)~ -(Sk + Sk+i)- 



(3) 



The factor j by the last term of Eq. (f3]l avoids a double counting 
of the Zeeman's term for the nodal Ising spins, which is equally 
split into two different cell Hamiltonians including one and the 
same nodal Ising spin. The partition function of the mixed spin- 
(1,1/2) Ising diamond chain can be written in this form 



N N 
lS k ) k=\ m ma lS k ) k=l 



(4) 



where fi = l/(ksT), k^ is the Boltzmann's constant, T is the 
absolute temperature and the symbol E(s t } marks a summation 
over all possible spin configurations of the nodal Ising spins. 
After performing the latter two summations over spin states 
of two interstitial Ising spins fi k ,i and fi k ,2 one gets the effec- 
tive Boltzmann's factor, which can be subsequently replaced 



with a simpler equivalent expression through the generalized 
decoration-iteration transformation 111 U Il2[ 1 1 3 1 



Zk ■ = ^ ^ exp(-/??4) = 2 exp 



— (Jk + Ok+l) 



P 



(J x +J 3 )(S k + S k+1 )-fih 



x jexpl — — - 1 cosh 
+ expl — — J cosh 
= A exp fiRS k S M + ^— (s 2 k + S 2 M ) + (3QS 2 k S 

+ pL(SkS 2 k+l +S 2 k S M ) + / ^(Sk+S k+ i) 



- (Ji - J3) (Sk - Sk+i) 



2 

A + l 



(5) 



As usual, the transformation parameters A, R, D, Q, L, and 
ho are given by a self-consistency condition of the decoration- 
iteration transformation (0, which requires that this mapping 
transformation must be valid independently of the spin states of 
two nodal Ising spins S k and S k +\ involved therein 01 llll2lll3l1 . 
By substituting all nine available spin states of two nodal Ising 
spins Sk and S k+i into the transformation formula (0 one 
merely gets six different expressions for the effective Boltz- 
mann's factor Zs k js M = Zk(S k , S k +i) 

Zo,o = 2 exp ) cosh (fih) + 2 exp (—^ j , 



Z±i,o = 2 exp I + y j|exp \ cosh 



f (h ~ h) 



cosh 



P 



^(J 1+ J 3 +2h) 



+ exp (-^) ( 

Zi,i=2 exp 03/ ! ){exp + ex p(~x) cosh [j3( J l+ J } -h)]^ , 

Zi,-x=2 exp j cosh [fi (/1-/3)] +2 exp cosh (fih) , 

Z_!,_i= 2 exp (-/?/*) {exp (^) 

+ exp (-^) cosh \fi(Ji+J 3 +h)]\, (6) 

which unambiguously determine so far unspecified mapping 
parameters through the relations 



A = Zo,0, /3R=\ln 
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(7) 



If the decoration-iteration transformation (|5) with the mapping 
parameters obeying the relations ((6)-(|7} is now substituted into 
Eq. ([8]) one in turn obtains a simple mapping correspondence 



Z(fi, Ji,J 2< h,h)= A n Zbeg(P, R,D,Q, L, h ), 



(8) 
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which relates the partition function X. of the mixed spin- 
(1,1/2) Ising diamond chain to the partition function 2^beg of 
the corresponding spin-1 Blume-Emery-Griffiths (BEG) chain 
Hill \m given by the effective Hamiltonian 

N N N 

k=l 



k=\ 



N 



k=\ 
N 



L kS l +1 + SfS k+i) - h ^S k- 



(9) 



k=\ 



k=\ 



It is quite evident from the effective Hamiltonian (0 that the 
mapping parameters R, D, Q, L, and ho represent the effective 
bilinear interaction, the single-ion anisotropy, the biquadratic 
interaction, the two-spin third-order interaction and the mag- 
netic field of the corresponding spin-1 BEG chain. 

The exact solution for the partition function of the general- 
ized spin-1 BEG chain can easily be found by means of the 
transfer-matrix approach 11151 Il6l Il7ll . Let us therefore merely 
recall the basic steps of this well-known procedure. The par- 
tition function of the spin-1 BEG chain can be first factorized 
into the following product 



.BEG - 

S\ So S n k=l 



(10) 



where the expression T(S k ,S k+ {) is defined as 

J3D 



T(S k ,S k+ i) = exp 



pRS k S k+1 + P — (S 2 k + S 2 k+l ) + PQS 2 k S 
+l3L{s k S 2 k+l + S 2 k S k+l ) + ^ (S k +S M ) 
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The relevant expression given by Eq. (fTTT i can be considered as 
the usual transfer matrix 



T(S k ,S k+ i) = 



( T(l,l) T(1,0) T(l,-1) N 
T(0, 1) T(0,0) T(0,-1) 
T(-1,1)T(-1,0) T(-l,-l), 



whereas the sequential summation over spin states of the spin-1 
BEG chain will correspond to a multiplication of the relevant 
transfer matrices. As a result, the partition function can easily 
be calculated with the help of the respective eigenvalues Ai of 
the transfer matrix with regard to 



Zbeg= J] T w (Si„Si) = TrT w = 2/if 



(12) 



Si=±l,0 



For completeness, let us quote the final expressions for the three 
transfer-matrix eigenvalues 



Ai = r + 2 sgn(g) -\fp cos 
where 



+ (l'-l)y 



(13) 



sgn(<?) 



-1, q <0 
k 1, q>0 



r = - [1 + 2 expOS^ + j3D + j3Q) cosh (J3h + 2/3L)] , 

1,1 , 
p = _( r - l) 2 + - exp(2/3/? + 2/3D + 2j3Q) sinh 2 (J3h Q + 2y3L) 

1 2 
+ - exp(-2/5R + 2/5D + 2J3Q) + - exp(J3D) cosh(f3h ) , 

q = r 3 - exp(J3R+2{3D+{3Q) cosh (2J3L) + exp(-J3R+2J3D+J3Q) 
+rexp(J3D) cosh (fih Q ) + (1-r) exp(2/3D + 20 Q) sinh(2yS/?) 
-rexp(/3R + (3D + PQ) cosh(ph a + 2/?L), 

1 ( 4p T:r q 2 \ 
0=-arctan — — . (14) 

3 { q ) 

In the thermodynamic limit — > oo, the partition function as 
well as the associated free energy per site of the spin-1 BEG 
chain is simply given by the largest transfer-matrix eigenvalue 

^max = max{A u A 2 ,A 3 } 



fi 



BEG 



-o lim 77 ln ^BEG = -4ln/l n 

P iv^cx, P 



(15) 



Our rigorous calculation for the partition function of the 
mixed spin-( 1,1/2) Ising diamond chain is thus formally com- 
pleted, since it is now sufficient to substitute the exact result 
( [Tol l for the partition function of the corresponding spin-1 BEG 
chain into the relevant mapping relation © between both parti- 
tion functions. As a result, the reduced free energy of the mixed 
spin-( 1,1/2) Ising diamond chain per unit cell reads 



/ = /beg - j InA = --(lnA + hU max ), 



(16) 



while the total magnetization normalized with respect to its sat- 
uration value readily follows from the relation 



m 

m. 



Idf 
2dh 



d In A <91n/l max 
+ 



d(ph) d(J3h) 



(17) 



3. Results and discussions 

Now, let us proceed to a discussion of the most interesting 
results obtained for the mixed spin-( 1,1/2) Ising diamond chain 
with the antiferromagnetic interactions (J\,J<2,Ji > 0), whose 
magnetic behaviour may be fundamentally affected by the ge- 
ometric frustration triggered by the competing diagonal inter- 
action J2 present in the diamond-like units. It is quite evident 
that the magnetic properties of the mixed spin-( 1,1/2) Ising di- 
amond chain given by the Hamiltonian ([D remain invariant un- 
der the transformation (JuJi,) — > {J^,J[), which allows us to 
consider /; > /? without loss of generality. For simplicity, 
we will pass to the dimensionless parameters JtlJ\, J3/J1, and 
h/J\ by normalizing all the interaction parameters with respect 
to the exchange constant J\, which will hereafter serve as the 
energy unit. While the former interaction parameter /2//1 mea- 
sures a degree of the geometric frustration inherent in the inves- 
tigated mixed-spin diamond chain, the latter interaction param- 
eter /3/./1 € [0, 1] characterizes a degree of the asymmetry of 
both Ising couplings along sides of the diamond units. 
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Figure 2: The ground-state phase diagram in the J2/J1 - h/J\ plane for: (a) the symmetric diamond chain with J-i/Ji = 1; (b) the asymmetric diamond chain with 
/3//1 < 1. The ground-state phase boundaries are explicitly quoted along the relevant lines in Fig. [2]). 



First, let us establish the ground-state phase diagram for the 
symmetric as well as the asymmetric version of the mixed spin- 
(1,1/2) Ising diamond chain. Depending on the interplay be- 
tween the interaction parameters Ji/Ji, J3/J1, and h/J\, one 
finds in total four different ground states to be further referred to 
as the antiferromagnetic state (AF), the nodal antiferromagnetic 
state (NAF), the unsaturated paramagnetic state (UPA), and 
the saturated paramagnetic state (SPA). The relevant ground 
states are unambiguously given by the following spin config- 
urations quoted along with the respective ground-state energies 
per primitive unit cell 



|AF>: 



n \Sk = -1) k 1 = 2) k 2 = 3) 

u\s k = i> ki = 4) K.2 = -5) 



h 



£af = —J\ + — - J3, 



I NAF) = 



n \s k = (-D") \i k A = HP) |mw = 
n \s k = (-D* +1 ) ki = k 2 = HH 



->NAF : 



T J 2 T 
-Jl ~ ~ + h, 



|UPA> = f]|S t - 1) 
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i> 



> 



N 1\ 1\ 

|SPA> = ]~||Sfc=l> ^,1 = 2) ^2=2) 

^2 

Sspa = 7i + — + ^3 - 2/z. 



(18) 



Apparently, the two-fold degenerate AF ground state corre- 
sponds to the antiferromagnetic ordering, at which the nodal 
Ising spins are aligned in opposite to the interstitial Ising spins. 
Even though the nodal spins have a twice as large magnetic 



moment as the interstitial spins, the total magnetization com- 
pletely cancels out on behalf of a twice as large number of 
the interstitial spins0 Another two spin configurations inherent 
to the NAF ground state involve the antiferromagnetic align- 
ment of the interstitial spins from the same diamond unit, as 
well as, the antiferromagnetic alignment of the nodal spins 
from the nearest-neighbour diamond units. Under this condi- 
tion, one observes an interesting doubling of the magnetic unit 
cell when comparing it with the primitive diamond unit cell of 
the mixed-spin Ising diamond chain. However, the most in- 
triguing spin alignment can be found in the macroscopically 
degenerate UPA ground state, where the nodal spins are fully 
polarized by the external magnetic field and the interstitial spins 
from the same diamond unit occupy one out of two equiprob- 
able antiferromagnetic states = 1/2)1/^2 = -1/2) and 
IjUjj = -1/2)|^ 2 = 1/2). The macroscopic degeneracy of the 
UPA ground state is accordingly proportional to a total num- 
ber of the interstitial spin pairs, which is also reflected in the 
respective value of the residual entropy *S res = M:Bln2. At 
sufficiently high magnetic fields, one finally detects the trivial 
SPA ground state with a full alignment of all nodal as well as 
interstitial spins into the external magnetic field. 

The relevant ground-state phase diagram including all the 
available ground states is displayed in Fig. [2] for the particular 
case of the symmetric diamond chain (Fig. |2^), as well as, the 
more general case of the asymmetric diamond chain (Fig. |2j)). 
Let us at first comment the ground-state phase diagram of the 
symmetric diamond spin chain. It is quite obvious from Fig.|2^ 
that the interaction term J2/J1, which is responsible for a geo- 
metric spin frustration, enhances a stability region of the UPA 
ground state, whereas there are two different scenarios of the 
magnetization process depending on whether J2/J1 < 4 or 



3 Note that both sublattice magnetizations completely cancel out just when 
assuming the equal g-factors for two different magnetic entities as originally 
imposed in the Hamiltonian 0. In the real magnetic substances, one should 
however expect at least some small difference between the relevant g-factors, 
which should result in a non-zero total magnetization and the ferrimagnetic 
rather than the antiferromagnetic spin arrangement. 
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Figure 3: (a) The total magnetization normalized with respect to its saturation value as a function of the magnetic field for the asymmetric diamond chain with 
J2/J1 = 4, /3//1 = 0.5 at a few different temperatures. Note that this choice of the parameters coincides with that one reported in Fig.4b of Ref. (TcJ; (b) 
Temperature variations of the total magnetization for the asymmetric diamond chain with J2 jJ\ = 4, J3/J1 = 0.5 at several values of the magnetic field. 



J2/J1 ^ 4. In the former case, the AF spin alignment is be- 
ing the respective ground state at low enough fields, the UPA 
ground state develops in a range of moderate fields and finally, 
the SPA ground state is stabilized at sufficiently high fields. On 
the other hand, the AF spin arrangement does not exist in the 
ground state for the latter particular case with fa/Ji > 4. Un- 
der this condition, the macroscopically degenerate UPA spin ar- 
rangement always represent the respective ground state below 
the saturation field, which determines the field-induced tran- 
sition towards the fully aligned SPA ground state. As far as 
the ground-state phase diagram of the more general asymmet- 
ric diamond chain (Fig.[2]3) is concerned, one observes the same 
general trends in the relevant ground-state phase diagram with 
exception of a presence of the additional possible ground state 
that corresponds to the NAF spin arrangement. It can be easily 
understood from Fig. |2j3 that the existence of the NAF ground 
state is restricted to the parameter region, where J 2 /J\ > 4/3 / J\ 
and simultaneously h/Ji < 1 — /3//1 . The former inequality im- 
plies that the NAF ground state replaces the AF one in a rather 
wide region of the parameter space when considering the highly 
asymmetric diamond chain with /3//1 << 1, which is also the 
condition of a sufficiently high persistence of the NAF ground 
state with respect to the external magnetic field according to the 
latter inequality. 

To provide an independent verification of the established 
ground-state phase diagram, let us examine in detail the magne- 
tization process of the mixed spin-( 1,1/2) Ising diamond chain 
by investigating the field dependence of the total magnetization 
at different temperatures. For illustration, Fig. depicts the to- 
tal magnetization normalized with respect to its saturation value 
in dependence on the reduced magnetic field for the fixed val- 
ues of the interaction parameters J2/J1 - 4 and J3/J1 - 0.5. It 
is noteworthy that the zero-temperature magnetization curve, 
which is in accordance with the ground-state phase diagram 
shown in Fig. [2J3, is plotted in Fig. [3^ by a thin solid line 
along with the finite-temperature magnetization curves calcu- 
lated with the help of Eq. ( fTTl i. In agreement with the ground- 
state phase diagram (Fig.|2j)), the zero-temperature magnetiza- 



tion curve exhibits two abrupt magnetization jumps reflecting 
two consecutive field-induced transitions between NAF-UPA 
and UPA-SPA at the relevant transition fields h c i/J\ = 0.5 
and h C 2/J\ = 3.5, respectively. It can be clearly seen from 
Fig. [3^ that the low-temperature magnetization curves closely 
follow the displayed zero-temperature magnetization curve (see 
for instance the magnetization curve for k B J '/J\ - 0.01), which 
proves a correctness of the established ground-state phase dia- 
gram and disproves an existence of any further ground state that 
would be reflected in some additional magnetization plateau. 
Apparently, the observed magnetization jumps and magneti- 
zation plateaus at zero and one-half of the saturation magne- 
tization are just gradually smudged upon increasing tempera- 
ture. To summarize, the low-temperature magnetization curves 
of the mixed spin-( 1,1/2) Ising diamond chain exhibit at most 
two different magnetization plateaus before reaching the satu- 
ration magnetization, which coincide with a presence of the AF 
or NAF ground state with zero total magnetization and of the 
UPA ground state with the total magnetization equal to a half 
of the saturation magnetization. 

Next, let us focus on typical temperature dependences of the 
total magnetization as depicted in Fig.|3j) for the same set of the 
interaction parameters J2IJ1 = 4 and /3/./1 = 0.5 at several val- 
ues of the external magnetic field. If the magnetic field is lower 
than the first critical field h < h c \, then, the total magnetiza- 
tion exhibits a round temperature-induced maximum as a func- 
tion of the temperature when falling towards zero at sufficiently 
low and high temperatures. Contrary to this, the total mag- 
netization starts from the one-half of the saturation value for 
the mediate magnetic fields h e (h c \,h C 2). The magnetization 
then either shows a monotonous temperature-induced decline 
for the magnetic fields slightly above the first transition field 
(i.e. h 7z h c \), or the non-monotonous temperature dependence 
with a round maximum for the magnetic field slightly below the 
saturation field (i.e. h iS h C 2). Finally, the total magnetization 
always exhibits a more or less steep temperature-induced de- 
crease when starting from the saturation value for the magnetic 
fields stronger than the saturation field h > h c z- For complete- 
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ness, let us quote that the total magnetization normalized with 
respect to its saturation value starts from the non-trivial values 
0.34151 and 2/3 for two particular magnetic fields equal to the 
transition fields h c \ and h C 2, respectively. 

Last but not least, let us make a few comments on the mag- 
netization curves of the mixed spin-( 1,1/2) Ising diamond chain 
obtained by Zihua Xin et al. II 1 QTI by employing the Monte 
Carlo simulations. Note that Xin and co-workers have reported 
in Ref. [ 10] two additional intermediate magnetization plateaus 
at 0.283 and 0.426 of the saturation magnetization, which were 
ascribed to a presence of some metastable states in the low- 
temperature magnetization curves. It should be pointed out, 
however, that the standard Monte Carlo simulation based on the 
Metropolis algorithm has been employed in Ref. ifioll . which 
should only give the stable states in thermal equilibrium rather 
than metastable states. From this perspective, one has to refute 
both striking magnetization plateaus at 0.283 and 0.426 of the 
saturation magnetization reported in Ref. Ill 0B . since they ev- 
idently contradict the exact analytical results presented in this 
work. It is worthwhile to remark that Monte Carlo simulations 
were performed in Ref. ifioll at the unusually low temperature 
k^T/Ji = 0.0001, which might indicate extremely long relax- 
ation times needed for establishing the thermal equilibrium and 
hence, one may consider an insufficient number of Monte Carlo 
steps used for equilibration at a given very low temperature as 
the main reason for the above mentioned discrepancy. 

4. Conclusion 

In the present article, the ground state and magnetization pro- 
cess of the mixed spin-( 1,1/2) Ising diamond chain have been 
rigorously studied by combining the generalized decoration- 
iteration transformation with the transfer-matrix method. In 
particular, our attention was focused on possible magnetization 
scenarios leading to the intermediate magnetization plateaus 
and the overall nature of available ground states. It has been 
demonstrated that the ground-state phase diagram of the sym- 
metric diamond chain totally consist of three different ground 
states, while the ground-state phase diagram of the more gen- 
eral asymmetric diamond chain includes in total four different 
ground states. It has been evidenced that two different mag- 
netization scenarios may in principle occur for the symmetric 
as well as asymmetric mixed-spin diamond chain depending on 
a mutual interplay between the interaction parameters. It has 
been actually evidenced that the low-temperature magnetization 
curves of the symmetric and asymmetric diamond chains dis- 
play at most two different magnetization plateaus, which mani- 
fest a presence of the AF or NAF ground states with zero total 
magnetization and/or UPA ground state with the total magneti- 
zation equal to a half of the saturation magnetization. Owing to 
this fact, our exact analytical calculations refute recent Monte 
Carlo simulations by Zihua Xin et al. 11 Oil , which have pre- 
dicted two additional striking magnetization plateaus at 0.283 
and 0.426 of the saturation magnetization. 

From the methodological point of view, we have adapted in 
the present work the generalized decoration-iteration transfor- 
mation in order to establish a rigorous mapping correspondence 



of the investigated mixed spin-( 1,1/2) Ising diamond chain with 
the effective spin-1 BEG chain. To the best of our knowledge, 
this form of the generalized mapping transformation has been 
adapted so far just for a calculation of the zero-field proper- 
ties of the mixed spin-(l,3/2) Ising linear chain [18, 19]. The 
present work thus brings a rather simple recipe to greatly extend 
this rigorous mapping approach, which may be further sub- 
stantially generalized in order to investigate magnetic proper- 
ties of several exactly soluble mixed spin-( 1,5) Ising and Ising - 
Heisenberg diamond chains in a non-zero magnetic field. As 
a matter of fact, the present approach can be rather straight- 
forwardly extended in order to account for the quite general 
spin numbers of the interstitial spins, the more general Heisen- 
berg interaction between the interstitial spins, the single-ion 
anisotropy, the next-nearest-neighbour and/or biquadratic inter- 
action between the nodal spins and so on. 
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